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Resume 

Nous donnons une representation stochastique pour une classe generale d’equations 
d’Hamilton-Jacobi (HJ) visqueuses, convexes et super-nonlineaires, au moyen d’equa¬ 
tions differentielles stochastiques retrogrades (EDSR) avec contraintes sur la partie 
martingale. Nous comparons nos resultats avec la representation classique en termes 
d’EDSR (super)quadratiques, et montrons notamment que I’existence d’une solution 
de viscosite a I’equation visqueuse de HJ pent etre obtenue sous des conditions de croi- 
ssance plus generales, incluant des coefRcients et une donnee terminate non bornees. 

Abstract 

We provide a stochastic representation for a general class of viscous Hamilton-Jacobi 
(HJ) equations, which has convex and super linear nonlinearity in its gradient term, 
via a type of backward stochastic differential equation (BSDE) with constraint in the 
martingale part. We compare our result with the classical representation in terms of 
(super)quadratic BSDEs, and show in particular that existence of a viscosity solution 
to the viscous HJ equation can be obtained under more general growth assumptions on 
the coefficients, including both unbounded diffusion coefficient and terminal data. 
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1 Introduction 


Given T G (0, oo) and d G N\{0}, we consider the parabolic semilinear partial differential 
equation (PDE) of the form: 


— — {t,x) — ^ti'[aa'^{x)D‘^u{t,x)] — b{x) ■ Dxu{t,x) 
+ F(^x,u{t,x), Q'^{x)Dxu{t,x)) 

u{T, x) 


0, (t,x) G [0,r) X M"*, (1-1) 

g{x)^ X G W^. 


Here a'^ and stand for the transpose of a and g, DxU and D^u are the gradient and the 
Hessian of u, respectively, and • represents the inner product between two vectors in 
The functions b: —)• a, ^—)• and F: x M x —)• M are called coefficients 

of the equation, and g: —)• M gives the terminal condition. 

The focus of this paper is to provide a probabilistic representation for solutions to (|l.ll) . 
This is a classical problem when g = a and the standard approach is to study a backward 
stochastic differential equation (BSDE) with the generator F and a forward diffusion X, 
namely: 


I xi’" = x + /; + /; a(X^’^)dWr, (t, x) G [ 0 , T] x t < s < T, ^ 

\ = g(X^^)-ffF(X^’^,Yj’^,Z*’^)dr-ffz*’^dWr, 

where VE is a d-dimensional Brownian motion. The forward backward stochastic differential 
equation (II.2p is then formally connected to the PDE (|l.ll) by the nonlinear Feynman-Kac 
formula: Y^'^ = u{t,x). Assuming that F is superlinear (which includes the quadratic 
and superquadratic case) and convex in its last argument (in which case, the PDE (II.ip is 
referred to as the generalized viscous Hamilton-Jacobi equation), existence and uniqueness 
of a solution to (|1.2p is a quite complicated issue, extensively studied in the literature. This 
paper provides an alternative representation, which does not require the classical structural 
condition g = cr, nor any non-degeneracy condition on g and a. In particular, our result 
can be applied to fi.t-„.d.ec,na«„n...e., = 00 Comparison to the literature will be 
presented in Remark 1 1.1 1 and in four examples below. To state our main result, we introduce 
the following assumptions which will be in force throughout the paper. 


Standing Assumption 1.1 

(i) There exists a constant such that 

| 6 (x) — 6 ( 3 :^)|-(-||(T(a;) — (t(x^)II-(-|| 0 (a;) — ^>(x^)|| < Lb^f^^g\x — x'\, for all x,x' 

where ||A|| = denotes the Frobenius norm of any matrix A. 

1. When (T = 0, the filtration F at the beginning of Section [2] below is generated by the Brownian motion 
B only. 
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(ii) There exist constants Mg > 0 and pg G [0,1] such that 

||^(x)|| < Mg{l + \x\^‘^), /or a//X G 

(iii) There exists a constant Lp such that 

\F{x,y,z) — F{x,y’,z)\ < Lp\y — y'\, /or any y, y'G M anc? (x, z) G x 

(iv) The map z F(x,y,z) is convex. There exist constants mp^Mp > t), q > p > 1, 
pp > 0, and qp G [0, (1 — Pg)q/{q — 1)) if Pg < 1, or qp = 0 if pg = 1, such that 

-mp[l + \x\P^ -l\z\P) < F(x,0,z) < Mp[l + \x\^^ + \\z\'i), 
for any (x, z) G x 

(v) Define the convex conjugate (or Fenchel-Legendre transform) of F as 

/(x, a,y) = — inf [a ■ z + F{x, y, z)], for all (x, a, y) G x x R. (1-3) 
The map x f{x,a,y) is continuous. 

(vi) The function g is continuous. There exist constants mg, Mg,qg > 0, and pg G [0, (1 — 
Pg)q/{q - 1)) ifPe < I, or Pg = 0 if pg = 1, such that 

—mg(l + |x|^®) < g{x) < Mg(l + \x\'^^), for all xGR'^. 

Let us now comment the above assumptions and their connection with related literature. 

Remark 1.1 1. The condition g > p > 1 in Assumption 11.11 (iv) means that F has a 
super linear growth in z. We also allow F to have some polynomial growth in x, and we 
distinguish the growth coefficients pp and qp for the lower and upper bounds. Indeed, notice 
that no condition is required on pp, while we impose some upper bound on qp depending 
on q and the (sub)linear growth coefficient Pg in Assumption 11.11 (ii). Observe that this 
upper bound qp = (1 — Pg)q/{q — 1) is decreasing with pg and q, with a limiting value 
equal to infinity when q goes to 1, and equal to 1 — when q goes to infinity; meanwhile 
1 — shrinks to zero (i.e. F is upper-bounded in x) when p^ = 1 (i.e. g satisfies a linear 
growth condition). Similarly, the terminal function g is enabled to satisfy a polynomial 
growth condition with the same constraint only on the power pg of the lower bound. These 
one-sided growth constraints on F and g are important for applications (see Example 11.31 
below) and are also sharp (see Example 11.41 below). 

When q = 2, F has at most quadratic growth in z. A stochastic representation of 
o with p = (T is given by the quadratic BSDE, which has been studied extensively, see 
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[m El Ena [251 [31 El [u, amongst others. For example, in the Markovian case of |17) . F and 
g are assumed to be bounded in x, i.e. qp = PF = Qg = Pg = 0, moreover g = a satisfies 
a linear growth condition, i.e. pg = 1. This case is covered by our assumptions (actually, 
we only need to assume qp = Pg = 0 when pg = 1, but no condition is required on pp, qg). 
In the Markovian setting of El, m or |25| . Q = a is assumed to be a bounded function, 

1. e. Pg = 0. Notice that we do not assume any non-degeneracy condition on a so that the 
case of time dependent coefficients as in |12| or |25], can be embedded in our framework by 
extending the spatial variables from x to (t, x). When p > 2 in Assumption 11.11 (iv), F has 
super-quadratic growth in z. This case has been studied in giEaiiaiiiiissiES], which will 
be compared with our results in Example 11.21 below. Together with three other examples, 
we shall illustrate in further detail the scope of Assumption 1.1 and compare our conditions 
on qF,Pg to existing results from both analytic and probabilistic aspects. 

2. The case where z F{x, y, z) is concave can be deduced from the convex case. Indeed, 
set F{x,y,z) = —F{x,—y,—z). Then F is convex in z and our results apply to equation 
a with F in place of F and —g in place of g. 

3. Assumption 1 1.1 1 1 vl is satisfied under the following two sufficient conditions: 

— The map x i—)■ F{x,y,z) is continuous, uniformly with respect to y and z, i.e., for all 
X G and e > 0, there exists <5 = 6{x,e) > 0 such that, 

\F{x,y,z) — F{x',y,z)\ < e, for any |x — x'| < (5 and (y, 2 ) G M x 

— The map (x, z) 1 —)■ F{x, y, z) is continuous for any y G M, and there exists an optimizer 

z* = z*{x,a,y) for (|1.3p such that z* is continuous in x. 0 

Example 1.1 (Generalized deterministic KPZ eqnation) Consider F{z) = —\\z\^ 
for some constants A > 0 and g > 1. The equation with Q = a as the identity 

matrix, is referred to as the generalized deterministie KPZ equation with the q = 2 case 

introduced by Kardar, Parisi, and Zhang in connection with the study of growing sur¬ 
faces. In this case, F{z) = satisfies Assumption II. 1[ in particular, pg = 0 in (ii) and 
Pf = qp = m. (iv). This equation has been studied from the analytical point of view in, for 
example, SI, in, and [15]. In particular, concerning Assumption [LT] (iv), in m Theorem 
2.6] inax{pg,qg} < q/{q — 1) is assumed, while here we only assume: pg < q/{q — 1), but 
no restriction on qg. Moreover, when max{pg,qg} = q/{q — 1), an example whose solution 
explodes in finite time is presented in |15l Remark 4.6]. Together with the uniqueness re¬ 
sult in |15l Theorem 3.1], this example shows that global existence of solutions satisfying 
\u{t,x)\ < C{1 -|- is in general not expected when pg = q/{q— 1). 0 
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Example 1.2 (BSDEs with superquadratic growth) Motivated by the previous exam¬ 
ple, Delbaen, Hu, and Bao HU studied BSDEs whose generator has superquadratic growth 
in the “Z-component". A solution to this superquadratic BSDE provides a probabilistic 
representation for the solution of (HU) with g = a and q > 2. In particular, given a forward 
process 

dXs = b{s, Xs)ds -|- adWs, 


with a Brownian motion W, a continuous differentiable function b with bounded derivative, 
and a constant matrix u, consider the following BSDE 

rT rT 


Ys = giXr) - [ F{Zu)du - f Z^dW^, 
J s J s 


where g is bounded and continuous, F is nonnegative, convex, and satisfies F{0) =0 and 
= oo. Proposition 4.4 in HU presents a solution to the previous superquadratic 


lim' ■ — 


BSDE. This result has been extended in |25) and HU, where F can depend on x and y, 
without convexity assumption on z, and cr can be a deterministic function of time. In these 
cases, Pq = 0, HU Assumptions (B.l) and (TC.l)] implies that max{pF,qF} < q/{q — l) and 
max{pg,qg} < q/{q — 1). When a depends on x, only existence for small time is available, 
see |25l Proposition 3.1]. Compared to aforementioned works, our main result (see Theorem 
11.11 below) provides an alternative representation for a global solution of (II.Ih when F is 
convex in z, p, not necessarily equals to cr, could depend on x and have (sub)linear growth. 
Moreover no restrictions are imposed on pF and qg. In particular, when p = a is bounded, 
i.e. Pg = 0, our result assumes only qF,Pg < 1- In fact, the asymmetry between upper and 
lower bounds of F and g is rather natural from BSDE point of view. Consider the BSDE 


= giWr) - 


1 I 7 
2\^u 


^du — 


ZudWu, 


whose solution is explicitly given by = — logE[exp(—(^(H t’))]. The previous expectation 
is well defined when g is bounded from below by a sub-quadratic growth function, however 
no growth constraint on the upper bound of g is needed. This asymmetry in assumptions 
also appeared in HU recently. <0 


Example 1.3 (Utility maximization) Our framework allows us to incorporate the two 
following financial applications. 

(i) Portfolio optimization. Consider a factor model of a financial market with a risk free 
asset and risky assets S = {S^, • • • , £'”) with dynamics 

dS^ = S^rdt and dSt = (Yiag{St)[{rln,-\-g{Xt))dt + a{Xt) dBt], (1-4) 

where r E M, /r and a are measurable functions on W^, valued respectively on M” and 
such that aa'^ is invertible, diag(S') is a diagonal matrix with elements of S on the diagonal, 
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In is a n-dimensional vector with every entry 1, and i? is a n-dimensional Brownian motion. 
We denote by A = a’'(oa'^)“^/r, the so-called Sharpe ratio. In the dynamics (|1.4I) . the factor 
X is a d-dimensional process governed by 

dXt = ^{Xt)dt + a{Xt)dWt (1.5) 

where /3 and a are Lipschitz functions on valued respectively on and and 

W is a d-dimensional Brownian motion. The correlation between B and W is given by 
d{W, B)t = pdt for some p G 

An agent with power utility [/(w) = vB for 7 < 1, 7 7 ^ 0 invests in this market 
in a self-financing way in order to maximize her expected utility of terminal wealth at 
an investment horizon T. Let v{t, w, x) be the value function of the investor and define 
the reduced value function u via v{t,w,x) = {w"^ Then the following equation, 
satisfied by u, is of the same type as (HU with Q = a (see e.g. [201 Equation (2.14)]): 

{ du 

— -^ — ^tT:[aa'^D‘^u\—b-DxU + F{x,a'^Dxu) = 0 

u{T,.) = 0, 

where 

b{x) = I3{x) + j^a{x)p^ X{x), F{x,z) = —^zMz^ — h{x), 

M = diag(ld)-b = 1 r + \j^\\{x)\^. 

We assume that b is Lipschitz. Moreover, note that M is positive definite. Hence T(x, z) = 
— F{x,—z) satisfies Assumption 11.11 (iv) with p = q = 2, and whenever A satisfies the 
condition: 

—mi?(l-b |x|^^) < 7 |Ap < Mf{ 1 + \x\‘^^), 

for some nonnegative constants mp, Mp, pp, and qp < 2(1 — Pg). Hence, when 7 > 0, 
such condition holds whenever A satisfies a strict sub-linear growth condition; while for 7 
< 0 (the empirically relevant case in financial context), the previous condition holds once 
A satisfies a polynomial growth condition. In particular, the second scenario includes the 
case where a constant, and p{x) (thus A(x)) is affine in x, as in the original Kim-Omberg 
model where X is an Ornstein-Uhlenbeck process with a constant and /3(x) affine in x. 

(ii) Indifference pricing. We consider a financial model as in (|1.4D - (|1.5I1 . where the process X 
represents now the level of nontraded assets (e.g. volatility index, temperature), correlated 
with the traded assets of price 5, for which the Sharpe ratio A is assumed to be bounded. 
Given an European option written on the nontraded asset, with payoff g{Xp) at maturity 
T, and following the indifference pricing criterion (see e.g. |19|L we consider the problem 
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of an agent with exponential utility U{w) = —7 > 0, who invests in the traded assets 
S' up to T where he has to deliver the option g{XT)- Let v{t, w, x) be the value function of 
the agent, and define the reduced value function u via: v{t,w,x) = U{w — u{t,x)). Then, 
u satisfies an equation of type ( 11 . 11 ) with g = a: 


du 

— — —— b ■ DxU + F{x,a'^Dxu) = 0 

uiT ,.) = g, 


where 


b{x) = /3{x) + a{x)p^X{x), F{x,z) = -^z{diag{ld) - p'^p)z'’+ jr + 


i|A(x)l^ 


Since A is assumed to be bounded, F{x, z) = —F{x, —z) satisfies Assumption ll.il (iv) with 
p = q = 2, pf = = 0. Moreover, Assumptions 11.11 (ii) and (vi) enable us to consider 

unbounded diffusion coefficient a and unbounded payoff function g, for example when X is 
governed by the “shifted" CEV model with <t(x) = ao{cri +xY‘>, for some positive constants 
do, di, and pg G [ 0 , 1 ) (the introduction of the positive constant di ensures that d is a 
Lipschitz function), and g satisfies the growth condition (recall Remark ll.11 2.1: 


—XIg{l + \x\'^^) < g{x) < mg(l + |xl^®), 
for some nonnegative constants nig, Mg, Qg, and 0 < Pg < 2(1 —pg). 
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Example 1.4 (A g ^ a case and beyond) The following linear-quadratic control prob¬ 
lem was studied by Da Lio and Ley m pp. 75 ]. Given valued deterministic functions 
A,B,C, and D, consider a linear stochastic differential equation 

dXs = [A(s)A, + B{s)as\ds + [C(s)X, + D{s)]dWs, Xt = x. 

Here a is taken from A which is the set of M'^-valued predictable measure controls. Given 
another M'^^'^-valued deterministic function Q, a constant matrix S G and a positive 

constant R, consider the following linear-quadratic (LQ) problem 


By noting that 


V(t,x)=infE / [XlQ(s)Xs + R\as\‘^]ds + XTSXT 
»eA L 


— inf \(DxuyB{t)a + R\a\‘^] = ynlBityD^ul 
aGM'* 


the Hamilton-Jacobi-Bellman equation associated to this LQ problem is written as 

f - A{t)x ■ DxU - x^Q{t)x + ^\B{tyDxu\'^ = 

\ u{T,x) = 


7 


0, on [0, T) X M'^, 
x'^Sx, x G M'^, 





where a{t,x) = {C{t)x + D{t)){C{t)x -\- D{t)y . In this case, a{t,x) = C{t)x + D{t) / 
B{t) = Q{t,x). 

More generally, equation was studied in |10) via viscosity solution techniques. 

Comparing to Assumption I l.li |10| restricts to bounded g, i.e., Pg = 0, and iJiax{pF,qF} < 
q/{q — 1). When max{pg,qg} < q/{q — 1), the existence of a viscosity solution to (|l.ll) is 
established in [TOl Theorem 2.2]. 0 


Before presenting our main result, let us first present an equivalent formulation of (HH) in 
terms of an Hamilton-Jacobi-Bellman (HJB) equation. Recall the convex conjugate function 
/ in ()1.3I) . By convexity of F in z, we then have the following duality relationship 


F{x,y,z) = — inf \a ■ z + f{x,a,y)], for all (x, y, z) € x M x ’ 
aeK'^ 

Therefore equation (ll.ip can be rewritten as the following HJB equation: 

— — {t,x)— inf \C°'u{t,x) + f(x,a,u(t,x))] = 0, (t,x)€[0,T)x 

ot 


u{T,x) = g{x), 


X G 


where 


1 


C°‘u{t,x) = {b{x) + Q{x)a) ■ Dxu{t,x) + -tT:[aa'’{x)D‘^u{t,x)]. 


( 1 . 6 ) 


(1.7) 


It is standard to relate equation (II.6p . at least formally, to the following optimal stochas¬ 
tic control problem of a recursive type: 

rT 


u{t,x) = infE 


as,u{s, Xl^-^^))ds + y(X^"’“) 


( 1 . 8 ) 


where the infimum is taken over all M'^-valued predictable processes a and the controlled 
diffusion evolves according to the equation 

= x + +Q{Xl'^^^)ar)dr + t<s<T, (1.9) 

where IT is a d-dimensional Brownian motion. 

Rather than studying the control problem (II.8p directly, following |16| . we introduce 
a randomized control formulation. For every {t,x,a) G [0,T] x x we consider the 
following forward system of stochastic differential equations: 


Xt,x,a = x + (6(Xr’“) + y(A:r’“) lr)dr + // aiXp^’^dWr, 


Is — a + Bg — Bt, 


( 1 . 10 ) 


where R is a d-dimensional Brownian motion, independent of W. In the next section, 
we shall check, under Assumption 11.11 (i) and (ii), that there exists a unique solution 





/*’“). System (11.101) is the randomized version of the controlled dynamics (|1.9I) . 
More precisely, the randomization procedure is performed by introducing the independent 
Brownian motion B, which is the natural choice when the control process a takes values 
in the entire space as in the present case. On the contrary, if the control process a is 
M-valued, for some compact subset A of then a natural randomization is carried out by 
means of an independent Poisson random measure /i on [0, oo) x A] see |16| . 

Now we introduce a stochastic representation for (II.6p and (|l.ll) . via the following BSDE 
with diffusion constraint. Given (t, x, a) G [0, T] x x consider 

n = r Yr)dr - {Kt - Ks) - f Z^dWr - VrdBr, (1.11) 

J s J s J s 

for all s G [t,T], together with the constraint 

Vs = 0, ds (8) dP-a.e. (1-12) 

The presence of the constraint (I1.12p forces the introduction of the nondecreasing process 
K. In Theorem 12.11 in Section [2] below, we construct and prove the existence of a unique 
maximal solution Yt,x,a^ (Il.lip - (I1.12I) . This allows us to present the 

main result of this paper, whose proof is presented in Section [31 

Theorem 1.1 For all {t,x) G [0,T] x and a, a' G 

_ Yt,x,a' 

Define 

u{t,x) := for all {t,x) G [0,r] x (1-13) 

Then there exists a constant C such that 

\u{t,x)\ < C{1 + for all {t,x) G [0,r] x (1.14) 

and u is a viscosity solution to equation (HI) ( or, equivalently, to (11.61) 1. 

Remark 1.2 Theorem 11.11 focuses on global solutions of (II.ip . i.e. T is not necessarily 
small. Local solutions have also been studied in the literature. When g is bounded, i.e., 
Pq = 0, and vaacicfpF, qp■,Pg-, <lg\ ^ dUd ~ 1)) [IHl Theorem 3.2] establishes a solution u 
satisfying 

|u(t, x)l < (7(1 + for some (7, (1-15) 

and equation dEH) on [T — T, T] for some r > 0. Similar local existence has also been 
obtained via probabilistic methods for (super)quadratic BSDEs in |121 Proposition 4.2] and 
|25l Proposition 3.1]. For global existence, the restriction on pg in Assumption 11.11 Ivil is 
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actually sharp. When Pg = 0 and Pg = q/{q — i), |10l Exemple 3.3] presents a deterministic 
control problem (i.e. cr = 0) whose value function is the only possible viscosity solution of 
(HH) in the class of functions satisfying (jl.lSp . However, under some parameter specification, 
the value function blows up when time is less than T — r for some r > 0. Therefore, similar 
to Example II. 1[ when pg = (1 — Pg)q/{q — 1), global viscosity solutions of (II.ip satisfying 
(ll.lSp do not exist in general, and the growth constraint in Assumption 11.11 (vi) cannot be 
improved even when dLH) is a first order equation (i.e., a = 0). 

Remark 1.3 When g is bounded, max{pF, qr^Pg, qg} ^ q/{q — 1), a comparison result for 
sub(super)solutions to (|l.ll) in the class of functions satisfying the growth condition (|1.14p 
was obtained in cni Theorem 3.1]. As a result, u given by (11.131) is the unique solution to 
equation dEU in the class of functions satisfying (II.lip . In this case, u is also continuous, 
as a byproduct of the comparison result. 0 


Remark 1.4 When there exists a smooth solution u to the equation (|l.ll) (or equivalently 
to the PDE (11.61) 1. one can check directly the connection between PDE (II.ip or (]1.6I) and 
the BSDE with diffusion constraint (|1.10p - (ll.lll) - (|1.12p . Indeed, by applying Ito’s formula 
to u along the forward diffusion process governed by (ll.lOp . we see that Yg = u{s, 
t < s < T, satisfies the relation (|l.lll) with 

,du 


Kg = 


[V 


dt 


+ C}-^' u) (r, + /(X*’^’“, 1 ^’“, y,) 


dr 


Zg = y, = 0, t<s<T. 


Since u is a solution to the PDE (|1.6I) . the term inside the bracket of K is nonnegative, 
which implies that {Kg)t<s<T is a nondecreasing process starting from Kt = 0 , and therefore 
the quadruple {Y, Z,V, K) is solution to the BSDE with diffusion constraint (ll.lOp - dl.lip - 
(I1.12p . Actually, this holds true whenever u is a subsolution to the PDE (II. 6 p . and the fact 
that u is a solution to this PDE will imply that (Y, Z, V, K) is the maximal solution to the 
BSDE with diffusion constraint (ll.lUp - (ll.lip - (|1.12p . This intuition will be proved in the 
rest of the paper without assuming the existence of a smooth solution to dnj. 0 


2 BSDE with diffusion constraint 

Let us introduce our probabilistic setting. Consider a filtered probability space (D, Xt, F = 
(J-'<j)o<s<T)F)) where F is the standard augmentation of the filtration generated by two d- 
dimensional independent Brownian motions W = (Ws)q<s<t and B = (i?s)o<s<T- For 
0 < f < T, we also consider the “shifted" version (D,J^,F* = (J'j)t<s<'r,F), where 
(tyi)o<s<T and {Bg)o<g<T before are replaced by {Wg - Wt)t<s<T and {Bg - Bt)t<g<T, re¬ 
spectively. We denote by E* the conditional expectation under P given y* for 0 < f < s < T, 
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and see that lE^ coincides with the “ordinary" expectation E. On this “shifted" probability 
space, we introduce the following spaces of stochastic processes. 

— §^(f,T): the set of real-valued cadlag E*-adapted processes Y = {Ys)t<s<T satisfying 


|y|P := E 


sup jy* 

t<s<T 


< OO. 


!^(f, T): the set of valued F*-predictable processes Z = {Zs)t<s<T satisfying 

rT 


IZI 




:= E 




n 


< OO. 


T): the set of nondecreasing F*-predictable processes K = {Ks)t<s<T £ T) 
with Kt = 0. We have 


\K\ 


's2(t,T) 


= E[K^ 


This section focuses on the construction of a maximal solution to the forward backward 
SDE with diffusion constraint (|1.10p - (jl.lip - (ll.l2l) . We first check the existence of a unique 
solution to the randomized forward system (jl.lOp . and show some useful moment estimates. 
Throughout this section, Assumption 1 1.1 1 is in force. 


Lemma 2.1 For any [t, x, a) £ [0, T] x x there exists a unique (up to indistinguisha- 
bility) adapted and continuous process {(As’^’“,Is’“),f < s < T} solving (ll.lOp . Moreover, 
for every m > \, there exists a positive constant C, depending only on m, T, Mg, 

and pg such that 


Ei 


sup 

re[s,r] 


< c(l + |A:*’^’“r+ / E^|/*’“li-"^]dr ), if Pg<l, (2.1) 


and, if pg = l. 


Et sup 
-re[s,r] 


< C I + IW 


t,x,a\m 


+ 


El[\Ir'^\^^]dr)yEl[e^fs ( 2 . 2 ) 


for all 0 <t < s <T. 


Proof. For every {t,a) £ [0,r] x there exists a unique process = {lt’°')t<s<T 
satisfying the second equation in (jl.lOp . Since the coefficients b{x) + Q{x)a and (j{x) of the 
diffusion system (ll.lOp for {X,I) are locally Lipschitz in (x,a), it is well-known (see e.g. 
Exercise IX.2.10 in m) that for every {t,x,a) £ [0,T] x x there exists an adapted 
process such that, if e := inf{s > t: = oo} with inf0 = oo, 

then jg unique (up to indistinguishability) adapted and continuous process on 

[t,e) n [t,T] that satisfies the first equation in (jl.lOp on [t,e) n [t,T]. It then remains to 
prove that the explosion time e of satisfies: E(e = oo) = 1. 
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For simplicity of notation, we denote I = and X = Consider T„ = inf{s > 

t: |Xs| > n} A T and apply Ito’s formula to for m > 1. From the dynamics (11.101) of 

{X,I) and Young’s inequality, we see, under Assumption 1 1.1 1 fif and (ii), that there exists 
a constant C (which in the sequel may change from line to line) depending only on m, T, 
Mg, and pg, such that for all t < s < T„, 

sup \Xrr < c(l + \Xsr + {\Xrr + \Irr + \Xrr-^^P‘>\Ir\)dr 

?’e[s,T„] V Js 


+ sup 


\Xur-^aiXu)dWuy (2.3) 


Case pg < 1. By applying Young’s inequality to taking conditional 

expectation on both sides of (|2.3p . and using the Burkholder-Davis-Gundy inequality, we 
get 


E! 


sup \Xr 
r&{s,Tn] 


< c(l + \Xr+ 


rp 

+ [ Es[ sup {\Xurii,^Tr,]{u))]dr 


which shows, by Gronwall’s lemma. 


E* sup |Y,r < c l + |Y,r+/ E‘[|/,|i-^.l[,, 7 ,„](r)]dr 

^rG[s,T„] -^ \ Js 

< c(^l + \Xsr +J\l[\lr\^]dry 

Since Tn e AT as n goes to infinity, from Fatou’s lemma we obtain 

pT 


E( 


sup \Xr 
rG[s,eAT] 


< C l + |Y,r+ / Ei[\Ir\^-^‘^]dr). 


(2.4) 


In particular, taking s = t, we have E(e = oo) = 1. Then, from (12.41) we deduce the required 
estimate (EH). 

Case Pg = 1. Take the conditional expectation in (12.3p with respect to the cr-algebra 

Ql := J-\ V cT(Ir) t < r < T), and observe that {Wg — Wt) t<s<T is a Brownian motion 

with respect to {Q\)t<s<T since W and I are independent. Therefore, we can still use 

t. G 

the Burkholder-Davis-Gundy inequality, and obtain (we denote by E^’^ the conditional 
expectation under P given Q\) 



sup 

■re[s,r„] 


X im 
y I 


< 


\IrrK'^ [hs,uir)]dr 

+ [ sup (|X„ri[,^T„](«))](! +l-^r-Ddr-Y 

Js uG[s,r] J 
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which gives, by Gronwall’s lemma and noting that T„ < T, 


E' 


t,G 


sup \Xr 

r£[s,Tn] 


< c(^l + \Xsr + \Irrdr^exp(^C (1 + |/^|)dr 

< C(^l + |X,r + \Irrdr^ d^dr_ 


Afterwards, by taking conditional expectation with respect to using the Cauchy-Schwarz 
inequality, and the fact that yja + h < y/a + y/h for any a, 6 > 0 , we get 


E* sup \Xrr < C{l + \Xsr + \ lKl[\Ir\‘^^]dr)JEl[e^fs dddry 
^re[s,r„] V Js ) 

Recalling that e f\T as n goes to inhnity, from Fatou’s lemma we obtain 


E*[|R|2™]dr WE^e^/^ (2.5) 


E‘ sup \Xrr < C l + |X,r + W 

^rG[s,eAr] \ V ' 

where the second conditional expectation on the right-hand side is finite, since / is a Brow¬ 
nian motion. Therefore, (|2.5p with s = t implies P(e = oo) = 1. Finally, from (12.5p we get 
the required estimate ( 12 . 2 p . □ 

Next, it is straight forward to check that Assumption 11.11 translate to the following 
properties on the generator / of the BSDE (jl.lip . 

Lemma 2.2 The map a f{x,y,a) is convex and satisfies 

\f{x,a,y)-f{x,a,y')\ < Lply-y'l, (2.6) 

-Mf{1 + \x\'^p - < f{x,a,0) < mpil + \x\p^ +-^\a\P'), (2.7) 

q Mp 

for all x € y,y' € M, and a £ Mfi. Here p' = p/{p — 1) and q' = q/{q — 1) are the 
conjugate exponents of, respectively, p and q in Assumvtion \l. 1[ 


We may then define the notion of maximal solution to the BSDE with diffusion constraint 

(frTT])-(fm. 

Definition 2.1 For every {t,x,a) £ [0,r] X X we say that 
iF*’^’“) £ S^(t, T) X ]Hl2(t, T) X ]Hl2(t, r) X ]K^(t, T) is a maximal solution to (ll.lll) - (|1.12p if it 
satisfies (Il.lip - (ll.l2p . and for any other solution £ §>'^{t,T) x 

e2(t,r) X H2(t,r) x K‘^{t,T) we have for any s £ [t,T]. 

Such a maximal solution is constructed using a penalization approach as in m Theorem 
2.1]. However, rather than employing an independent Poisson random measure as in |16| . 
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our randomization here is carried out by means of an independent Brownian motion in 
(ll.lOp . Let us consider, for every {t,x,a) S [0,T] x x and n G N, the following 
penalized BSDE: 


Y^i _ 


+ f - {K^ - K2) - j ZyWr - [ V^dBr 

J s J s J s 

+ r V^dr - f ZyWr - f K^dBr, ( 2 . 8 ) 

J s J s J s 


for all s G [t,T], where 

K = 

and the generator/„(x, a, y, u) = f{x,a,y) — n\v\. Notice bv (|2.6I) that this generator f„ is 
Lipschitz in (y, v), so that from standard result due to |21) . we know that for every {t, x, a) G 
[0, r] X X and n G N, there exists a unique solution yn,t,x,a'^ ^ 

S2(t,r) X X El2(t,r) to BSDE (ITKD . 


Moreover, we have the following comparison results. 


Lemma 2.3 For every {t, x, a) G [0, T] x x the following statements hold: 

(i) The sequence nonincreasing, i.e., > y?i+i,i,a:,a^ n G N. 

(ii) For any solution ^ ^ ]g 2 ^ jj 2 ^ 

to (fLTTTl - lfmi . we have > yb^.a^ n G N. 

Proof. Since /„ > fn+i, the first statement follows from a direct application of the compar¬ 
ison theorem for BSDEs. For the second statement, note that 1*’°', Yr’^’^, Vr’^’^) = 

/(Xr’^’“,due to Vr’^’^" = 0, and is a nondecreasing process. Then the 

second statement readily follows from [231 Theorem 1.3]. □ 


The aim is to obtain a uniform bound on (||T"'’*’^’“||§ 2 (j 71 ))^, which, together with the 
monotonicity property stated in Lemma D.dl iL allows to construct y*>^>“ as the limit of 
(y”’*’^’“)„. In contrast to [161 Lemma 3.1], where the compactness of the space of control 
actions A (which does not hold in the present case, since A = M'^) is exploited to prove 
the §^-bound, we utilize a dual (or randomized) representation of To this end, 

we introduce some additional notations. Let denote the set of i3„-valued predictable 
processes u, where Bn is the closed ball in with radius n and centered at the origin, and 
define T> = VJnBn- For t G [0, T] and 1 / G P, define the probability measure on 
via 


"dP 


= 8 


VgdBs ) = exp 


(/ 


i/g dBg 2 


ji^sl^ds 


where the Doleans-Dades stochastic exponential on the right-hand side is a martingale due 
to the definition of T). In the sequel, we denote by the condition expectation under P*^ 
given Jf, for any s G [t,T], and u gV. 
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Remark 2.1 Notice by Girsanov’s theorem that W remains a Brownian motion under 
for any v ^T). Then, by the same argument as in the derivation of (12.11) in Lemma [2.11 we 
see that when < 1, for every m > 1, there exists a positive constant C, depending only 
on m, T, ^ind pg such that 


sup < 

rGlSjT] 


c(l +|X*’^’“r+ / 


(2.9) 


for all {t,x,a) S [0,T] x x s G [LT], and v ^V. 


0 


The next result provides a dual representation of the solution to the penalized BSDE. 


Proposition 2.1 For every {t,x,a) G [0, T] x x and n G N, it holds that 

rT 


Yn,t,x,a ^ ess inf E*/ 




J s 


( 2 . 10 ) 


for all t < s <T, where the predictable process y"': fl X [t, T] —)> R is given by 

jt^a y"n,t,x,a^ _ jt,a 


7" = 




Y'n,t,x,a 

In particular, | 7 "’| is bounded uniformly by the constant Lp in (j2.6j) . 

Proof. Applying Ito’s formula to e-^s and integrating between s and T, we 

obtain 

T 

Yn,t,x,a ^ -n / 

J S J S 

rT rT 

_ / - / e/-’’ 

J 5 ./ S 

Take G Pru then add and subtract the term '^^'^'^Vf''^’^’°'Urdr, so that we obtain 

T 

yn,t,x,a ^ + J 0)dr 

_ e^s 7llrf«(^|yn,t,x,a| ^ Vff'^^^^^Vr)dr (2.11) 

fT 

- / e^s'fudnYn,t,x,a^^,^_ / ^lidnyn,t,x,a^QU^ 


where I's ds + B and W are both P'^-Brownian motions by Girsanov’s theorem. 

Note that both local martingales in (|2.1ip are in fact P'^-martingales. Indeed, from Bayes 
formula and the Gauchy-Schwarz inequality, we have 



< 


L'^ W \\Z\\ 


15 













where is bounded in L^(ri, J^,P) since v is bounded. Combining the previ¬ 

ous estimate with the Burkholder-Davis-Gundy inequality, we obtain that ZrdWr is of 
class (D) (see for instance Definition IV.1.6 in |24| 1 under P*^, hence it is a P'^-martingale 
(Proposition IV.1.7 in |24|1. The same argument can be applied to the other local martin¬ 
gale. Projecting both sides of (|2.11l) on J-”*, and using the inequality n|u| -|- u • a > 0, for 
any u £ M and a £ Bn, we obtain 

Yn,t,x,a Y ess inf 




J: 0)dr + T'- 


du 




To prove the reverse inequality, denote V”’’* the f-th component of £qj. ^ _ 

1,... ,d. Set i/*’* = — for r £ [t,T] and i = 1,... ,d. Then u* £ Vn, 
and = 0 for any r £ [f, T]. It then follows from (12.111) that 

■T 

e/: Jt,a^ ^ 'l^du 




from which the claim follows. 


□ 


Combining the above dual representation together with the moment estimate for X in 
Remark O we obtain the following uniform estimate. 


Corollary 2.1 There exists a positive constant C, independent of n, such that 

_j_ IJf*i2;,a|ijirVpg^ ^n,t,3;,a 

^ ^ PF'^’ig 

< 1 , 


Pp'Vqg f 

C{1 + ), 


c(^l + |V^’"^’“|Pr’V<?9 + |4-“|PFV5gVp'^gC|/i’“| ifp^ = i^ 


( 2 . 12 ) 


for all {t,x,a) £ [0,T] x x s £ [t,T], and n £ N. In particular, £ §^(t,T) 

and 

{ £Fy9£„ / 

C(1 + \:,\PFV,FVPgV<ig + |a| ^fp^ < 1, 

(5(1 -F Ixl^r’VgFVpgVgg |^|pp.VggVp'^gC|a|^ if ’ 

for some positive constant C, independent of n. 


Proof. The upper bound in (I2.12h . From the dual representation formula (I2.10jl . we 
have (recall that, when n = 0, P^ coincides with P) 

rT 


Yn,t,x,a Y 


e/s" l^duf^xt,x,a^ jt,a^ ^ 


t,X,a'y 

T ) 


Then, recalling the fact that \^'^\ is bounded by Lp, exploiting the polynomial growth 
condition of / in (12.7j) and of g in Assumption 11.11 (v), we get 


Y 


n.t.x.a 


< (C l + E 


sup + sup \II^\P' 

LrG[s,T] ?-g[s,T] 


(2.14) 
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Together with the estimate (|2.1I) . and the standard estimate E* [sup^gj^ < (7(1 + 

|jt,a|m) £qj^. g^j^y ^ 1, this sliows the upper bound in (|2.12p when < 1. 

When = 1, by using the estimate (|2.2I) . together with the fact that ab < (a^ + 6^)/2 
for any a, 6 G M, we obtain from (I2.14p : 


•yn,t,x,a ^ 
^ s — 


C 1 + 


|j^t,a:,a|pj?Vgg _|_ ^jt,a^ppVqgVp''^ [e*" 


(2.15) 


For the expectation on the right-hand side, we have 


_ gCT|/a’ Ie |^gC'^®“Po<t<T l-(^t I j ... E , (2.1 


where the last equality follows from the independence of B^,..., B’^, the components of the 
d-dimensional Brownian motion B. Now, notice that 


E ^^P0<t<T l-(^t I j . . . E l-Sf*! j 

... Esupg<t<j.(B^)+j ^/^E™po<t<T(-Bt)~j (2.17) 

Using the property that, for every j = l,...,d, {—B^)t>Q is still a Brownian motion, 
we see that the stochastic processes {{Bj)^)t>o = (max(i?^,0))j>o and {{Bj)~)t>o = 
(max(—0))t>o have the same law. As a consequence, the random variables supo<t<T(7?t )"*" 
and supQ<j<2’(7?i)~ have the same distribution. Moreover, the distribution of supQ<j< 2 ’(i?j )^ 
(or, equivalently, of supo<t<'r(-Bt)~) is independent of j = 1,.. ., d. Therefore, supo<4<'r(7?t )^ 
and supQ<j< 7 ’(i?j)“ have the same distribution as supQ<^<ji(i?j )^. Then, (|2.17p becomes 

E l-Sf I j ... E l-Sf*! j < E 


As supo<t<T7?t > 0, P-a.s., (since Bq = 0, P-a.s.), it follows that, P-a.s., supo<t<T(-Bt )'*' = 
supo<t< 7 ’• In conclusion, (|2.16p becomes 

E*[gt^/Jh^’“M-] < gCr|4-|E[g2CTsupo<,<^Bi]y (2.18) 


Using the reflection principle (Proposition III.3.7 in |24|i. and in particular that supQ<j< 2 ’ B^ 
has the same law as we obtain from (|2.18l) 

E*[gt^/Jd*’“M-] < gCr|/^“|E[g2CT|Bi|]'^ < ^gCT|/i'“| < ^gC|/i’“|_ (2.19) 
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Combining the previous two estimates (|2.15p and (12.191) . we confirm the upper bound in 
(I2.12P for pg = 1 case. 

The lower bound in (|2.12p . From formula (|2.10p . the polynomial growth condition of 
/ and g, we find that jg greater than or equal to the following quantity 


ess inf 




In the case where pg = 1, we have Qf = Pg = 0, and since Iq”"! is bounded by Lf, we then 
obtain the lower bound: > —2(MfT + mg)e^^'^. When pg < 1, we use the estimate 

()2.9I) and the fact that | 7 "’| is bounded by L^?, to deduce that there exists a positive constant 
C, depending only on qF,Pg,T,mf,q', Lh^cr,e, Lf, Mg and pg, such that 

Yn,t,x,a Y —C{l + 


+ ess inf E*’'^ 




M 


l-q' 


^—LpT\Yq' _ 


Recalling that qF,Pg < (1 — Pg)q', we see that the function h{x) = 

’IF _ 

C\x\ — C\x\ , X G M“, is bounded from below on M“ by a constant independent of n. 
Therefore, we obtain the lower bound. 

Estimate (12.131) . When < 1, it follows directly from bounds (|2.12p . Lemma [2T] 
with s = t, the standard estimate EPsup^gjj j^] < (7(1 + lal™) for any m > 1. 

When Pg = 1, (I2.13P follows from the Cauchy-Schwarz inequality together with I] < 

( 7 gC|a| latter inequality follows from EPe^l^®’ I] < = e^l“lE[e^l^®-*l] < 

QC\a\^'^^C afterwards we reason as in (12.161) 1. □ 

The previous uniform norm estimate implies the following uniform norm estimate on 

^n,t,x,a j^n,t,x,a'^ 

Corollary 2.2 There exists a positive constant C, independent of n, such that 


ryn,t,x,a\\ I \\'\7"n,t,x,a\\ i ii r^n.t.x.a 

^ ’ IL2(gr) + l|L ” ’ + 11-^ 


(t,T) 


ls 2 (t.T) 


< 


Pp'^’lg V, I 

(7(1 + _g 1^1 1-p, Vp 


( 2 . 20 ) 


ifPe < 1 

(7(1 + |a;|PFVgfVp9Vgg _g \a\PFVqg\/p'yC\a\^ 
for all {t,x,a) G [0,T] x X n G N. 

Proof. By proceeding along the same arguments as in the proof of m Lemma 2.3], we 
have: 

||^n,£,3i,a||2 _|_ ||'^n,t,3i,a||2 _j_ ||^n,£,3i,a||2 
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< C E 




+ sup || y ’"’*’^’“||2 
nGN 


Then, by exploiting the polynomial growth condition of / in (12.7p combined again with 
Lemma 12.11 for s = t, the standard estimate E[sup^g[j < (7(1 + |a|”*) for any 

m > 1, and I] < (7e^l“l when pg = 1, and using the estimate (I2.13h for we 

obtain the required uniform estimate. □ 


Now the previous uniform norm estimates allow us to take limit as n —>■ oo. 


Theorem 2.1 Let Assumvtion \l.l\ hold. For every {t,x,a) G [0, T] x x there exists 
a unique maximal solution (yi:a;,a^ g §2(t,T) x M‘^{t,T) x ]HI^(t,r) X 

]K^(t,T) to (ll.lll) - (|1.12p . such that, as n ^ oo, \ yh^’'^- j^n,t,x,a converges 

to in L^(n,J7j,P), for all s G [t,T]; and {^z^d,x,a^yn,t,x,a^^ weakly converges to 

^Zt,x,a^yt,x,a'^ in (t,T) X (t,T). Moreover 

{ PF^^g \ / f 

C{1 + + 1^1 i-p, Vp ^ ^ 

(7(1 + \x\PP'^iP^P9'^is + |a|Pf’V''9Vp')e^l“l, ifpg = 1 

where C is the same constant as in (|2.13p . 

Proof. The proof follows the same passages of the proof for [161 Theorem 2.1], with 
some small modifications due to the fact that here we adopted a Brownian (rather than 
Poisson) randomization. For this reason, we just outline main steps of the proof. For 
every {t,x,a) G [0,T] X x M'^, from Lemma [T3f il and estimate (I2.13p . it follows that 
the sequence converges decreasingly to some F^adapted process satisfying 

the bound (|2.2ip . Next, the monotonic limit theorem (see |23[ Theorem 2.4]) implies the 
weak convergence stated in the theorem and that the limit (y*>^>“^ ^ 

§^(t,T) X ]H[^(t, r) X ]HI^(t,r) X ]K?{t,T) satisfies BSDE (II.IIP In order to prove the con¬ 
straint ()1.12p . define the functional F{V) = Efor V G ]HI^(t,T). Notice that 
jg^n,t,x,a _ ^y^yn,t,x,ay Hence estimate (12.201) implies sup„ E|iLy < oo, therefore 

lim„_>.oo y(y”> 0 a:,a^ _ g g^j^^e F is lower semicontinuous with respect to the weak topol¬ 
ogy of M‘^{t,T) {F is convex and strongly continuous), it then implies F{y^’^’°‘) = 0 and 
()1.12p holds. Finally, the maximality of (yi:a;,a^ follows from a direct 

application of Lemma 12.3l fiil and the uniqueness follows from [231 Proposition 1.6]. □ 


3 Feynman-Kac representation formula 

The present section is devoted to the proof of Theorem 11.11 Assumption 11.11 is in 
force throughout this section. Let us hrst recall the definition of (discontinuous) viscosity 
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solution to equation (HSl) (or, equivalently, to m)- Define the lower semicontinuous 
(Isc) envelope u* and upper semicontinuous (use) envelope u* of a locally bounded function 
u: [0, T) X —>• M as follows: 

n*(t,a:)= liminf u{t',x') and u*{t,x) = limsup u{t\x') 

(t ,x it',x')^{t,x) 

t'<T 

for all (t, x) G [0, T] x 

Definition 3.1 

(i) We say that an use {resp. Ise) function u: [0,T] x —>■ R is a viscosity subsolution 

{resp. supersolution) to (11.61) if 

u{T,x) < {resp. >) g{x), VxGR'^ 

and 

— -;^{t,x)— inf \C°‘ip{t,x) + f{x,a,u{t,x))] < {resp. >) 0 
Ot aeK'^ 

for every {t,x) G [0,T) x R"^ and every ip G T] x R'^) satisfying 

{u — p){t,x) = max {u — (p) (resp. min (rt —w)). 

[O.TlxR-^ ^ [ 0 ,T]xK<i ’ 

(ii) We say that a locally bounded function u: [0,T) X R'^ —)■ R is a (discontinuous) 
viscosity solution to (USD whenever u* is a viscosity subsolution and is a viscosity 
supersolution to ([ISD- 


3.1 Penalized BSDE and corresponding semilinear parabolic PDE 

Let us define, for every n G N, functions Vn, v: [0, T] x R'^ X R"^ —^ R via 

Vn{t,x,a) = and v{t,x,a) = for (t, x, a) G [0, T] x R'^ x R*^. 

Notice that Vn v pointwise as n goes to infinity. Moreover, it follows from (I2.13P and 
()2.2ip that 


Vn{t,x,a)\ + \v{t,x,a)\ < 


PF'^igx, / 

C{l + + |o| ^ ), 

C{l + |3;|PFVgFVpgVg9 _|_ |^|PFVcjgVp'^ gClal ^ 


(3.1) 

if Pg = 1 


for all (t, X, a) G [0, T] X R'^ x R'^ and n G N. For each n G N, let us consider the following 
semilinear parabolic PDE 
' dv 

X, a) - ^AaX„(t, X, o) - C°'Vn{t, X, a) 

'• — f{x.,a,Vn{t,x,a)) + n\DaVn{t,x,a)\ = 0, (t,x, o) G [0,T) x R'^ x R'^, (3-2) 

Vn{T,x,a) = g{x), (x, a) G R'^ x R'^, 
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where is the Laplace operator with respect to a and is given in (|1.7p . Then, we have 
the following result. 


Proposition 3.1 The function Vn is a continuous viscosity solution to (|3.2p . i.e., Vn = 
(vn)* = {vn)* and properties in Definition \3. 1\ are satisfied where the equation is replaced by 

( 1 ^ . 


Proof. When pg < 1, so that Vn satisfies a polynomial growth condition, the result is well- 
known, see e.g. [221 Theorem 4.3]. When = 1, a change of variable is needed in order to 
deal with the exponential growth of in a. More precisely, define the map fi: M —?■ M as 
follows 


/3{x) 


- 1 , 
< Fix), 


X > log 2, 

— log 2 < X < log 2, 
X < — log 2, 


where P{x) = Ax^ + Bx^ + Cx (for some real constants A, B, C) is a polynomial of fifth 
degree which realizes a smooth paste of the two other branches of the map (3, so that 
fi G C2(M) (since P is an odd function, in order to determine A, B, C it is enough to realize 
a C^-paste with the upper branch of /3) and /3 is increasing, hence invertible. We denote by 
a: M —)■ M the inverse map of (3, which is given by 


a{y) = < 


log(l + y), 

P-Hv), 

^-log(l -y). 


y > 1, 

y < -!• 


(3.3) 


For a = (ai,..., Orf), 6 = {bi,...,b^) G we denote a{b) = {a{bi),..., a{bd)) G M'^, 
/3(a) = (/3(ai ),..., (3{ad)) G and define 


Wn{t,x,b) = Vn{t,x,a{b)). 

From m, with = 1, we obtain that there exists a positive constant c such that 

\wn{t,x,b)\ < + + (3.4) 

for all {t, x, b) G [0, T] x x and n G N. Notice that 

DaVn{t,x,a) = 




dbi 

d‘^Wn{t,x,fi{a)) , 


dbd 




2 = 1 


2=1 
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Then, Vn is a continuous viscosity solution to (|3.2D if and only if Wn is a continuous viscosity 
solution to the following equation: 

X, b) - ^A^Wn{t, X, b) - X, b) 

^ —f{x,a{b),Wn{t,x,b))+n\D^Wn{t,x,b)\ = 0, (t, x, 6) € [0, T) x x (3-5) 

Wn{T,x,b) = g{x), (x, 6) G X 


where 


n/3 u\ fdwn{t,x,b) , 

D^Wn[t,x,b) = ( - — -/? (a(6i)),... 




d‘^Wn{t,X,b) , 


dbd 

d 


A^Wnit,x,b) = + Y 


dwn{t,x,b) „ 


dbi 




2=1 ^ 2=1 
Proceeding as in m Theorem 4.3] we can prove that Wn is a continuous viscosity solution 
to the previous equation, so the claim follows. □ 


3.2 The invariance of the function v with respect to the variable a 

We distinguish between the two cases pg < 1 and pg = 1. Let us consider the two 
following hrst-order PDEs: 

\Dav{t,x,a)\ = 0, (t,X,a) G [0,T) X X (3-6) 

\Dbw{t,x,b)\ = 0, (t, X, 6 ) G [0, T) X X M'^. (3.7) 

Here 'w{t,x,b) = v{t,x,a{b)), where a{b) comes from (13.3p . 

Lemma 3.1 When pg < 1 (resp. pg = 1) the function v [resp. w) is a (discontinuous) 
viscosity solution to (13.6p (resp. (13.7p ). 

Proof. We firstly suppose that pg < 1, so that Vn and v satisfy a polynomial growth 
condition as reported in (13.ip . Since the supersolution property clearly holds, let us focus 
on the subsolution property. To this end, since x is a pointwise limit of the decreasing 
sequence (vn)n, it is upper semi-continuous and v = v*. Take (t,x,a) G [0,T) x x 
and ip G C'^’^([0,T] x x satisfying (v — p)(t,x,a) = maxp — p). From 

the polynomial growth property (|3.1I) . we can assume without loss of generality (up to a 
polynomial perturbation of p for large values of x and a) that the previous maximum is 
strict. Moreover, using standard techniques in viscosity solutions (for similar arguments, 
see for instance [2] or Section 6 in ID), we can deduce the existence of a bounded sequence 
(tn,Xn,an)n S [0,T) X X R'^ such that 

(vn - p)(tn,Xn,an) = max (vn - p) and 

[0,T]xR‘ixK‘^ 
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{t, X, a, v{t, X, a)) 


as n —>■ oo. 


{tn ) ) On j Vn (tn , Xn , fly;,) ) 


The viscosity subsolution property of Vn then implies that 

1 / d(p 1 

I Xn 1 Un) I ^ ~ f i^ni Xn j ®n) T Xni ®n) T ^ ^(tn: Xn: ®n) 

T f (^Xn: Ojn: Vn{tn: Xn: CLn}^ 

Sending n to infinity, and using the continuity of the coefficients b, a, and /, we deduce the 
claim. 

Suppose now that = 1. Reasoning as in the case < 1 and using the fact that 
P' > 0, we can prove that, for every i = 1 ,..., d, rc is a viscosity solution to the first-order 
PDE: 

= 0, (t, X, 6) G [0, T) X X 


dw{t, X, b) 


dbi 

As a consequence, we deduce that re is a viscosity solution to (13.71) . □ 

We can now state the following result on the independence of v with respect to a. 


Proposition 3.2 The function v satisfies 

v{t,x,a) = v{t,x,a'), for all {t,x) £ [0,T] and a,a' 

Proof. Let us firstly consider the case p^, < 1. Then, this result can be proved proceeding 
as in m Lemma 3.4 and Proposition 3.2], with only two small modifications. First, in 
|16| the equation —\Dav{t,x,a)\ = 0 is considered rather than (|3.6p (this is due to the fact 
that the function v therein is the increasing limit of {vn)n: since a Hamilton-Jacobi-Bellman 
equation with a “sup” operator, instead of “inf”, is studied there). However, the proofs of 
|16l Lemma 3.4 and Proposition 3.2] are not affected (apart from obvious modifications by 
the presence of the minus sign). Second, in |16) the variable a belongs to an open, bounded, 
and connected subset A of some Euclidean space rather than to the entire space 
as here. To link to |16| . let us consider the open ball Br C for any r > 0. Lemma l3.ll 
implies that u is a viscosity solution to equation (|3.6I) on [0,T) x x Br, for every r > 0. 
At this point we can apply [161 Lemma 3.4 and Proposition 3.2] to deduce that v satisfies 

v{t,x,a) = v{t,x,a'), for all {t,x) G [0,T] x and a, a' G Br- 

As a result, the claim follows since r is arbitrarily chosen. 

Suppose now that p^ = 1. Then, proceeding along the same lines as in the case p^ < 1, 
we conclude that w does not depend on b. Since v{t,x,a) = w{t,x, 13(a)), this implies that 
V does not depend on a. □ 
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3.3 Prom the BSDE with diffusion constraint to the viscous HJ equation 


Following Proposition 13.2[ we define the function u, as in (I1.13p . which satisfies (due to 
()3.1I) with a = 0) 

\u{t,x)\ < C{1 + for all {t,x) G [0,r] x 

Let us now prove that u is a viscosity solution to equation (II. 6 p (equivalently, to equation 
O), which, together with Proposition 13.21 concludes the proof of Theorem 11.11 

Proof of Theorem II.11 In order to prove that u is a viscosity solution to equation (|1.6I) . 
we begin by noticing that, as a direct consequence of Proposition 13.11 Vn is a viscosity 
solution to (j3.2l) on [0,T] x X Br, for any r > 0. Therefore, as a decreasing limit of {vn)n, 
V, hence u, is a viscosity solution to the following equation: 

du 

— -;^{t,x)— inf \C‘^u{t,x) + f{x,a,u{t,x))] = 0, (t, x) G [0, T) x 

< “ 6 ®- (3.8) 

u{T,x) = g{x), x G 

The proof of this result can be done proceeding as in US Section 3.4], with some small 
modifications, due to the Poisson, rather than Brownian, randomization used in m, which 
affects in part the form of the penalized PDF (13.21) . From the definition of viscosity solution 
to (|3.8p we have, for any {t,x) G [0,T) x M'^, r > 0, and any ip G C'^’^([0,T] X R*^), 

— -^{t,x)— inf [C°'p{t,x) + f{x,a,u{t,x))] < (resp. >) 0 

ut a^Br 

whenever (u — (/?)(t, x) = maxp 7 nj^]^d(ri — lyj) (resp. minjQ 2 ’]xRd(^‘~‘P))- Noticing that the 
choice of the test function p is independent of r, we deduce 

inf [C°‘p{t,x) + f{x, a, u{t,x))] < (resp. >) 0 , 
at ogR'^ 

which implies that u is a viscosity solution to equation (II. 6 p . Finally, the terminal condition 
is also satisfied thanks to the terminal condition in (13.21) and the fact that u is a monotone 
limit of {vn)n- n 
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